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An analysis of the influence of the Landau-Migdal-Pomeranchuk (LPM) effect on the development 
of air showers initiated by astroparticles is presented. The theory of Migdal is studied and compared 
with other theoretical methods, particularly the Blankenbecler and Drell approach. By means of 
realistic computer simulations and using algorithms that emulate Migdal's theory, including also 
0^ ' the so-called dielectric suppression, we study the behavior of the relevant observables in the case of 

0^ ' ultra-high energy primaries. We find that the LPM effect can significantly modify the development 

of high energy electromagnetic showers in certain cases. 
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'nJ" ! I. INTRODUCTION 

(N ■ 

^ ■ The study of atmospheric showers initiated by high energy astroparticles plays a central role in contemporary cosmic 
' ray physics The most important component, in number of particles, of the air shower is by far the electromagnetic 
one. Additionally, many shower observables that give direct information about the properties of the primary depend 
CO strongly on its behavior. 

Electron bremsstrahlung and pair production are the dominant processes in the electromagnetic shower at very 
high energies. The standard description of these processes, i.e., the Bethe-Heitler equation can be incomplete at 
qI^ ^ very high energies due to some effects that drastically reduce the corresponding cross sections ^ . Those mechanisms 
play a relevant role in the development of air showers because they can lengthen them, and consequently move the 
position of the shower maximum deeper into the atmosphere. 

In the present work we study two of the several suppression processes that can affect the high energy electromagnetic 
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O ' interactions ||] , namely, the Landau-Pomeranchuk-Migdal (LPM) effect and the dielectric suppression . The 
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first of these effects is due to the multiple scattering while the second one is due to the interaction of the bremsstrahlung 



^ ■ photons with the atomic electrons in the medium, through forward Compton scattering. 



The LPM effect was studied semi-classically by Landau and Pomeranchuk [|j and the first quantum mechanical 
approach was given by Migdal JH . Recently, the experimental confirmation of the LPM effect at the Stanford Linear 



Accelerator Center (SLAC) [g|,[7| has originated new theoretical works [g-lS] and several analysis of its consequences, 
particularly in cosmic ray physics p^ . 

The main purpose of our work is the study of the modifications in the development of the air showers that take 
place when the LPM effect and the dielectric suppression are taken into account. In order to do that we first study 
the different theoretical approaches of the LPM effect. In particular, we compare the Migdal formulation |^ with 
other approaches that were developed more recently [pP-[l0|] . We also analyze the modifications of the bremsstrahlung 
and pair production probabilities in the case of our interest, that is, when the medium is a "layer" of air of infinite 
thickness and variable density. Our study is, therefore, complementary to the recently published review by Klein 
1^ which primarily treats the LPM effect in the case of solid targets of finite size, like the ones used in the already 
mentioned SLAC experiment. 

To complete our study we have developed new LPM/dielectric suppression procedures, and we have installed them 
in the AIRES air shower simulation system ||l5|. The AIRES code has then been used as a realistic air shower 
simulator to generate the data used to make our analysis of the influence of the LPM effect on high energy showers. 

It is worthwhile to mention that AIRES represents a set of programs to simulate atmospheric air showers and to 
manage all the associated output data. Most of the physical algorithms of the AIRES system are based on the realistic 
procedures of the well-known MOCCA program 0]. AIRES provides some additional features, for example: The 
Earth's curvature is taken into account allowing safe operation of all zenith angles; the simulation programs can be 
linked to different alternative hadronic collision models like QGSJET or SIBYLL |18 1^; etc. 



The version of AIRES used in this work is capable of processing the following particles: photons, electrons, positrons, 
muons, pions, kaons, eta mesons, nucleons and antinucleons, and nuclei up to Z = 26. Among all the physical processes 
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that may undergo the shower particles, the most important from the probabihstic point of view are included in the 
simulating engine. Those processes can be classified as follows |]l5| : (i) Electrodynamical processes: Pair production 
and electron-positron annihilation, bremsstrahlung (electrons and positrons), knock-on electrons {S rays), Compton 
and photoelectric effects, LPM effect and dielectric suppression (discussed in this work), (ii) Unstable particle decays, 
pions and muons, for instance, (iii) Hadronic processes: Inelastic collisions hadron-nucleus and photon- nucleus, 
generally simulated using the external packages]^ (QGSJET or SIBYLL). Photonuclcar reactions, elastic and inelastic, 
(iv) Propagation of particles: All particles -even the short lived ones- are properly tracked taking always into account 
the corresponding mean free paths for the processes they can undergo.^ Continuous energy looses (ionization), multiple 
Coulomb scattering and geomagnetic deflection are taken into account for charged particles. 

This paper is organized as follows: In section |l| we introduce the different approaches to the LPM effect and 
compare the Migdal formulation with other ones. In section |ll^ we explain the practical implementation of the LPM 
effect and the dielectric suppression into the AIRES program. In section [V we analyze the results of the air showers 
simulations performed with such code. Finally, in section ^ we place our conclusions and comments. 



II. THEORY OF THE LPM EFFECT 



The so-called LPM effect takes place when an ultra-relativistic particle emits low energy bremsstrahlung photons 
passing through matter. In this case, fewer photons are emitted than predicted for the case of isolated atoms H. In a 
similar way, the cross section for electron-positron pair production is reduced in the case of high energy gamma rays 
passing through matter. 

This effect was first predicted by Landau and Pomeranchuk Q some 40 years ago. They treated the classical 
radiation of a high-energy particle in a fluctuating random field inside an infinitely thick medium. Afterwards, Migdal 
1^ provided the corresponding quantum mechanical theory giving analytical expressions for the bremsstrahlung and 
pair production cross sections when matter is present. Migdal theory was developed for an infinite target thickness. 
More recently, an experiment performed at SLAC [Q measured the LPM effect finding that there is an acceptable 
agreement between the experimental data and Migdal theory which is presently considered the standard treatment. 

The experimental confirmation of the LPM effect triggered new theoretical works |]^-p^. In particular, Blanken- 
becler and Drell Blankenbecler 10 1 and Klein ||] have reanalyzed the problem, upon better approximations 
than in the earlier works, and they were able to consider effects of finite target thickness. We include in this paper a 
comparative analysis between different approaches. 

To make our paper self-consistent, we start with a qualitative description of the LPM effect explaining the main 
features of Migdal theory. We then compare the final Migdal results with those of Blankenbecler and Drell. 



A. Bremsstrahlung 



Let us consider first the case of bremsstrahlung where an electron or positron of energy E and mass m emits a 
photon of energy k in the vicinity of a nucleus of charge Z . 

Neglecting the photon emission angle and the scattering of the electron, the minimal longitudinal momentum 
transfer^ to a nucleus, q||, in the case when E ^ k and E ^ mc^, is given by ^ 

km?c^ 
^11 = 2E{E-ky 



^The external hadronic packages simulate nontrivially the physical processes involved in a hadron-nucleus inelastic collision. 
Among other features, both QGSJET and SIBYLL are capable of generating high energy gammas coming from the decay 
of baryonic resonances, even if the probability of such diffractive processes diminishes as long as the energy of the projectile 
particle is increased |E0|. 

^This is particularly important for the case of electromagnetic showers initiated by hadronic primaries. The electromagnetic 
component comes mainly from 7r° decays, and even if they are short lived particles, the probability of undergoing hadronic 
collisions cannot be neglected at ultra-high energies. 

^The minimal longitudinal momentum transfered to the nucleus occurs when the final electron and the photon follow in the 
same direction of motion than the initial electron. 
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The formation length can be defined according to the uncertainty principle: Ifo = h/q\\. Ifo is the distance over which 
the interaction amplitudes can add coherently. It is straightforward to show ^ that if other interactions are present 
while the electron is traversing this distance, then the resulting amplitude will, in general, be reduced. 

Let us consider the case when multiple scattering takes place during the radiation emission process. Using a small 
angle approximation, the longitudinal momentum transfer can be expressed as 



2E{E-k) 2c 



M5/2 



where 9ms/2 is the multiple scattering angle in half the formation length, that is {Es/E)y/lf/{2Xo), with Eg = 

mc^^^-Kja = 21.2 MeV (a ~ 137-^), and Xq = [^■qar'^Z'^ ln(184Z-i/3)] is the radiation length (rg = e^/mc^ and 
77 is the number of atoms per volume unit). The effect of multiple scattering becomes significant when the second 
term of equation is comparable in magnitude with the first one. This is the case when 

^ E + E'lpm ^ ^ 

where 

» = |, (4) 

and 

is a characteristic, material dependent, energy which gives the scale where the effect cannot be neglected. For a 
given energy E, the emission of photons with y < E/{E + -Elpm) will be affected by the interference due to multiple 
scattering; and therefore for E ^ -Elpm the effect extends to the entire photon spectrum < y <1. 
The characteristic energy Slpm can be expressed as 

^™ " Wnhrjr^Z^ ln(184Z-i/3) ' ^ ^ 

This equation makes evident the fact that i?LPM diminishes when the density of the medium increases. Therefore, 
for dilute media, the LPM effect will be appreciable only for energies much higher than the typical ones used in 
experiments with dense targets In fact, for air in normal conditions E'lpm — 2.2 x 10* GeV; one can compare 

this energy with the energy of the electron beam of 8 and 25 GeV used at SLAG . 

The formation length, when multiple scattering is the dominating process, can be put in the form 



^'-^'"^ E{E-ky 

To measure the strength of the effect it is convenient to introduce the suppression factor through 

^ _ da-LPu/dk 
doBH/dk 

where dasH/dk stands for the bremsstrahlung cross section given by the classical theory of Bethe and Heitler 
The fact that the cross sections are proportional to the formation length gives rise to 



For k small in comparison with E, The cross section found by Bethe and Heitler [E| is proportional to 1/fc. If 
multiple scattering is taken into account, this proportionality changes, due to equation (0), to l/\/k. 

In the Migdal theory of the LPM effect Q, the multiple scattering is treated as a diffusion mechanism. The 
average radiation per collision and the interference between the radiation from different collisions are then computed. 
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When collisions occur too close together, destructive interference reduces the radiation. The multiple scattering is 
treated using the Fokker- Planck technique to solve the corresponding Boltzmann transport equation. Migdal includes 
quantum effects, such as electron spin and photon polarization, but his calculations only apply for a target of infinite 
thickness. The resulting cross section for the bremsstrahlung process reads: 

da,,u iarUis) ^ ^[1 + (1 - In ( '-^] (10) 



dk 3k \ Zs 

where 



]J 8E{E - k)as)' ^ ' 

2 if s < si 

£.{s)={ 1 + lns/ln.si ifsi<s<l, (12) 
1 if s > 1 



(yil = ^1/3/184) and 



G(s) =48s2|j+Im*(^s + i+is^|, (13) 

$(s) = 12s^ jlm [*(s + is) + ^{s + l + is)] - ^| : (14) 

where represents the logarithmic derivative of the complex F function. The functions G(s) and $(§) are plotted in 
figure |l|. Both G{s) and <i>(s) belong to the interval [0, 1] for all s > and their limits for s going to zero (infinity) 
are (1). Therefore, when s ^ 1 the suppression is important (the cross section is reduced, S 1), while for s 1 
there is no suppression [S = 1) and the Migdal cross section reproduces the main terms of the Bethe-Heitler equation 
if complete screening is considered, namely 

It is easy to see, for example, that if the medium is air^ at the altitude h which corresponds to a vertical depth 
Xy = 1000 g/cm^ {Xy{h) = J^^ p{z)dz H, Xy{0) = 1000 g/cm2, X„(100 km) = 0) and for electron energies of 100 
TeV, the screened]^ Bethe-Heitler equation and the no suppression limit ( p^ of Migdal theory agree within the 2.5% 
relative error level. In concordance with equation ^ the region of photon spectrum where the LPM suppression is 
significant grows when the primary energy is enlarged, and affects virtually the complete range < y < 1 when the 
energy of the electron is 10^* eV (in this case E^lpm = 2.8 x lO" eV = 280 PeV). 

The plots in figure ^ illustrate the dependence of the LPM effect with the density of the medium. In this figure the 
probability 

is plotted versus y for several electron energies. The medium is air at two representative altitudes: (a) Vertical 
atmospheric depth X^ — 1000 g/cm^, which corresponds approximately to an altitude of 300 m.a.s.l {p = 1.19 



*For altitudes up to 90 km above sea level, the air is a mixture of 78.09% 20.95% O2 and 0.96% other gases [g2j which 
can be adequately modeled as an homogeneous substance with atomic charge and mass numbers Zctt — 7.3 and AcS — 2 x ZcB, 
respectively. 

^The screening effect of the outer electrons has been calculated by Bethe and Heitler on the basis of the Fermi- Thomas model 
of the atom. The influence of the screening on a radiation process is determined by the quantity 7 = 100{m/E)[y /{l — y)]Z~^^^ . 
For 'y ^ 1 screening can be practically neglected while one has complete screening for 7 = |2U2]J]. 
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kg/m^, E'lpm — 280 PeV); (b) Xy = 50 g/cm^, which corresponds approximately to 20,000 m.a.s.l (p = 78 g/vc?, 
i^LPM = 5.78 EeV). The plots in figure ||a show that an important suppression takes place at all the energies 
considered and becomes severe for the 10^" eV case. On the other hand, the LPM effect does not seriously affect the 
bremsstrahlung probabilities for 10^^ eV electrons at — 50 g/cw? as shown in figure ||b. In this case the critical 
energy where the effect becomes significant is placed around 10"'^^ eV since i?LPM is some 20 times larger than the 
corresponding value for figure ^a. 

When the effect of the polarization of the medium is taken into account, a significant alteration in the bremsstrahlung 
formula for soft photons appears. The interaction of a photon with the atomic electrons produces another kind of 
suppression of the cross section that is called dielectric suppression If we take into consideration that the 

dielectric constant of the medium is different from one, i.e., e = 1 — {Uujpf' / k"^ , where ujp is the well-known plasma 
frequency (w^ = Ai:Ze^rj/m] for air in normal conditions hojp = 0.72 eV), then the longitudinal momentum transferred 
to a nucleus changes from equation to 

'^ii^2£;(£;-fc)+^^- ^^^^ 

Consequently, when the second term in the last equation is comparable with the first, the dielectric suppression 
becomes important. This happens when the energy of the photon is much lower than a critical energy given by 

fccrit = ^iWcrit with 



EiE-k) E , , 

Wcrit = WpV - 2 1^ 

It is worth to notice that the correction due to the dielectric constant is certainly negligible for the propagation 
phenomenon. On the contrary, in our case of the emission process, the effect is measurable because here the scale 
fixing parameter is the critical frequency given by equation (|l8|), and not the plasma frequency. Notice also that this 
critical frequency is exactly the plasma frequency when measured in the electron rest frame. If we define now 

hujp 

2/diel = 2' (19) 



equations ( |17D and ( |18D tell us that the dielectric suppression takes place when y <^ 2/dioi- 

The Migdal approach takes into account this dielectric effect. Since dielectric suppression occurs for y ^ y^id ^ 1, 
the term in G{s) can be neglected. In that case, the cross section becomes 

da^PM 16«r,2e(4^^,^)^,j /184\ s^l+fy^V. (20) 



dk 3kS ' ' \zi 

Figure || illustrates the influence of the dielectric suppression on the bremsstrahlung cross section (20) where the 
probability for bremsstrahlung is plotted against y in the case of i5 = 10^* eV, and for an altitude of 300 m.a.s.l. The 
emission probability is suppressed for y <C ydici — 1-4 x 10~^. This corresponds to photon energies fc ^ 1 TeV. 



B. Pair production 



The pair production process can be treated similarly as bremsstrahlung. The corresponding expression for the cross 
section of this interaction, where a photon of energy k produces a pair e'^e~ of energies E and k — E respectively, can 
be obtained immediately from the bremsstrahlung formula [^^. In this case the cross section reads 

^'^^^"^y^"^'^ = '-^{GiD + 2[u + (1 - ur]m}Z^ In (^) (21) 



where 

E 

and 



kEu 



PM 



8E{k - E)^{Tj 



(22) 



(23) 
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Clearly the dielectric suppression on the cross section of the electron-positron pair production can always be neglected. 
This follows from the fact that the polarization of the medium influences only the soft photons whose energies are 
much smaller than that of the electrons . 

Figure H shows the normalized probability of pair production at atmospheric depths of 1000 g/cm^ (a) and 50 g/cm^ 
(b) for different energies of the primary photon. The production probabilities are progressively suppressed when the 
primary energy rises. From figure ^, it is also evident that symmetric processes {u ~ 0.5) are the most affected by 
the LPM suppression. 



C. Comparison with other approaches to the LPM effect 

We should remember that the Migdal approach does not include all the corrections that should in principle be 
included. In fact, the inelastic form factor that accounts for inelastic interactions with the atomic electrons, and the 
term that accounts for the Coulomb corrections because the interaction takes place with the electron in the Coulomb 
field of the nucleus are not taken into account. However, Migdal formulation proves to work sufficiently well 

when its results are compared with experimental data . 

Blankenbecler and Drell ]^-pX|] have computed, in an alternative way, the magnitude of the LPM suppression. They 
have used the eikonal formalism, standard in the study of scattering from an extended target. This approach leads to 
a physically clear quantum-mechanical treatment of multiple scattering and from that to the derivation of the LPM 
suppression of soft photon radiation from high energy electrons in matter. The approach includes as limiting cases 
the Bethc-Hcitlcr radiation from a charged particle scattering against an isolated atom, relevant for a thin target 
and, in the opposite limit, the LPM effect which suppresses the radiation in a thick target. Their result for thick 
targets reads 

k-jr= (24) 



dk nEh 



where £ is the target thickness and 



w [ ^ '^^^ — 1 + ^wx — 1) ^.^ ^ , 

Jo 1 + Awx - VI + ^wx 

sin a; , tt , , 

-^dx^- (26) 



where 



2(1-2/) 



(27) 



E fl-y\ e , , 

6-E/LPM \ y J Xo 

An additional corrective term due to the phase-amplitude correlations has been calculated in reference , and must 
be added to the probability of equation (p^). In the case of a finite thick target, this term can be expressed as follows 

,dP\ 3a(l-w)r f'"/2wbi + l \, 

'^Jr^ie^l [—z lj(-l)^eosa:d., (29) 



with bi = i/£fo and z = y/l + Awx. 

The transition to the LPM regime occurs for u; = 1 and the extreme LPM limit for w ^ 1 in concordance with the 
relation (||). 

Migdal theory contains a number of approximations that are not very transparent on physical grounds On 
the other hand, the Blankenbecler and Drell formulation is theoretically more robust in its principles, but does not 
include dielectric suppression. Migdal results are simpler to treat numerically and the differences between the two 
approaches are small in the case of infinitely thick media. This shows up clearly in figure ^ where both predictions 
and their relative differences are plotted for representative energies. 

To make these plots we have evaluated numerically the integrals of equations ( p5| ) and (|9|). Our results can be 
considered as a generalization of the analysis presented in reference |^] which corresponds to the limiting cases of 
negligible (u; ^ 1) or severe {w ^ 1) suppression. 
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In reference [Q, other approaches to the LPM effect plj-ITst are reviewed; and similarly as in the case of the 
Blankenbecler and Drell formulation their results are not in contradiction with Migdal theory. It is worth noticing 
that the main differences between such different approaches appear when considering finite thickness slabs, which is 
not the case of interest in our work, devoted to the study of the LPM effect in an infinite dilute medium like the 
Earth's atmosphere. 



III. PRACTICAL IMPLEMENTATION 



During the computer simulation of a ultra-high energy air shower, the processes of emission of bremsstrahlung 
photons and e~ pair production need to be calculated some millions of times. This requires that the algorithms 
used in such calculations should be fast enough to let the simulations complete in a moderate amount of time. The 
Migdal formulation for bremsstrahlung and pair production can be implemented fulfilling those requirements [p^ and, 
as discussed in the previous section, its results are in acceptable agreement with other alternative approaches. To fix 
ideas, we shall illustrate the practical implementation in the case of bremsstrahlung; the algorithm for pair production 
is completely similar. 

The probability of equation ( [l^ ) can be put in the form 



where 



dk 3fc 



(31) 



and 

F{E, k) = ^{s){y'G{s) + 2[1 + (1 - yfMs)} (32) 



s, ^(s), G{s) and <I>(s) are defined in equations ( |ll[^14| ). Therefore it is possible to simulate the bremsstrahlung 
processes in two steps (i) Rough approach u sing the probability (|T]), (ii) Correction to (i) using the rejection 

approval algorithm to give the exact distribution (|30|). The correction factor is the function F{E,k) adequately 
normalized. Step (i) relates with the normal, Bethe-Heitler bremsstrahlung and can be processed straightforwardly. 
To complete step (ii) it is necessary to evaluate i^LPM, s, ^(s), G{s) and <&(s). £^lpm can be evaluated directly 
from equation (|^) and requires the estimation of the local air density p — —dXy/dh. s and ^(s) can be conveniently 
calculated by means of repeated iterations of equations ( |ll|) and (O) star ting with the initial value ^(s) — 1. This 
iterative process normally converges in a single step. Equations (|l3| ) and (|lj) allow to obtain accurate estimations 
for G{s) and $(5) using the standard algorithms to evaluate the complex function |2^, but are not suitable for 
production procedures which must be fast. In this case it is better to represent G{s) and $(s) by means of rational 



functions: 



^i^) = . ''/'''"-^y 3 (34) 

^ ^ 1 + diS + d2S^ + C3S^ ^ ' 



where the coefficients are adjusted to minimize the maximum relative error for s > 0. Notice that equations ( |33 
and ( |34| ) give exact results for both s = and s — > cx) The fitted coefficients are listed in table 1; the maximum 
relative error for G(s) ($(s)) is 0.15% (0.12%). 

The dielectric suppression can be easily implemented as a scaling factor to s, as it clearly follows from equation 
(pO|). In AIRES this scaling factor correction is applied when ydici/y > 1/2, i.e., the dielectric suppression is neglected 
in all cases where j/dici/y < 1/2. ydici is evaluated at each case taking into account the altitude of the particle. 
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IV. SIMULATIONS 



We have investigated the modifications in the bremsstrahlung and pair production cross sections due to the LPM 
effect, for individual processes in different conditions. In this section we are going to analyze the influence of the effect 
in the development of air showers initiated by ultra-high energy astroparticles. 

There are two characteristics that must be taken into account when analyzing the shower development: 

1 . The atmosphere is inhomogeneous, the density of the air diminishes six orders of magnitude when the altitude 
goes from sea level to 100 km , and therefore the characteristic energy i?LPM of equation (||) varies accordingly. 
As a result, the suppression in the corresponding cross sections depends strongly on the altitude where they 
occur. In figure I the energy E'lpm is plotted against the vertical atmospheric depth, for Xy in the range 0.1 
g/cm^ (66 km above sea level) to 1000 g/cm^ (about sea level). To illustrate the meaning of this plot let us 
consider a ultra-high energy electromagnetic process with a primary energy of, say, E = 10^" eV = 100 EeV: 
If the primary particle is located at an altitude of 100 g/cm^ the process will be strongly suppressed whereas 
if it is located at < 1 g/cm^ the suppression is negligible. In a similar way, the dielectric suppression will 
depend on the altitude since j/dici depends on the air density. However, this parameter does not change with Xy 
so dramatically as i?LPM (compare figures |^ and 

2. The number of particles in the ultra-high energy showers is very large (about 10^^ for a primary energy of 10^*^ 
eV), and therefore the influence of the LPM effect and/or dielectric suppression on global observables may not 
be visible if the fraction of affected particles is not large enough. 

To estimate how the LPM effect can modify the shower development when the preceding features are taken into 
account, we have run a set of computer simulations using the AIRES program [T5| ] to determine what is the average 
fraction of electromagnetic particles (gammas, electrons and positrons) significantly affected by the LPM suppression. 
We have performed simulations using 3 x lO^'' eV gamma or proton primaries. The particles were injected at the top 
of the atmosphere (100 km above sea level) with vertical incidence. 

At various predetermined altitudes (observing levels), the fraction, />, of particles whose energies are greater than 
-^'LPm/4 were recorded for each simulated shower. This fraction is taken equal to 1 before the first interaction takes 
place.^ The data related to every one of the defined levels were then averaged. In figure I the average fractions, (/>), 
are plotted versus the vertical depth. The primaries are gammas. The results coming from simulations where the 
LPM effect is taken into account are plotted (full line) together with similar data obtained without evaluating such 
effect (dashed line) . We can see that in both cases the relative number of particles suffering the LPM effect diminishes 
when Xy grows, reaching almost zero for points that are always above 500 g/cm^. 



Of course, the plots represented in figure 
at a different altitude, the fraction of partic 



3 must be regarded just as qualitative indicators. Since each shower starts 
es capable of undergoing the LPM effect will vary accordingly: If a shower 
starts very high in the atmosphere, at Xy < 5 g/cm^ for example, it is likely that the LPM effect will not significantly 
affect its development since the first interactions will take place in a region where i?LPM is very high (see figure]^). 
On the other hand, a penetrating primary particle starting the shower at, say, Xy > 100 g/cm^ will surely show a 
significant modification in its longitudinal development due to the LPM suppression. 

It is worthwhile mentioning that the influence of the LPM effect on the shower development depends also on the 
inclination of the shower: For large zenith angles it is more probable that the showers will start at points where the 
density of the air is very low, and so a smaller suppression should be expected in this case. 

In the case of proton primaries, the electromagnetic shower starts after the first hadronic interactions (mainly after 
7r° decays) and so the energies involved are smaller than the primary energy. For this reason, (/>) vanishes faster 
than in the case of gamma primaries. In fact, this fraction is virtually zero for Xy > 100 g/cm^ (for simplicity we 
have not included any proton related data in figure 

The most evident signature of the impact of the LPM effect on the shower development is the shift on the position 
of the shower maximum, ATmaxJ^ This can be easily seen in figure || where the total number of charged particles 
for showers initiated by gamma rays is plotted against the vertical depth. These data were obtained by means of 



^In the case of 3 x 10^" eV gamma showers (LPM effect taken into account), for example, the mean value of the first interaction 
depth is approximately 63 g/cm^ with a standard deviation of 80 g/cm^. These figures become respectively 46 g/cm^ and 45 
g/cm^ when the showers are simulated without considering the LPM effect. 

'^The shower maximum, Xmax, is defined as the atmospheric depth of the point where the total number of charged particles 
is maximum. The number of charged particles at Xmax is noted as A'^max. 
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computer simulations using the program AIRES jist , and correspond to inclined showers with zenith angle 60 degrees; 
the primary particles are injected at the top of the atmosphere and the ground level is located at a vertical depth of 
1000 g/cm^. Two primary energies are considered, namely, 10^^ eV (figure ^) and 3 x 10^" eV (figure ^b). 

The severe suppression suffered by the first interactions that start the shower shows up clearly in figure ^ where 
the plot corresponding to the simulations performed taking into account the LPM effect shows significantly more 
penetrating showers that develop deeper in the atmosphere in comparison with the non-LPM ones. The position of 
-'^max is shifted in about 130 g/cm^ (about 260 g/cm^ measured along the shower axis) when the LPM procedures 
are switched on. The number of particles at the maximum, -/Vmax is also affected: It is about a 40% smaller in 
the LPM suppressed showers. This is mainly due to the fact that the electrons and positrons loss more energy 
in ionization processes before being able to generate bremsstrahlung photons, and therefore the shower reaches its 
maximum generating less particles than in the non-LPM case. 

For lower primary energies the influence of the LPM effect on the development of the shower is less significant, as 
shown in figure ^ a for the case of 10^^ eV gamma showers. Accordingly with the results of our simulations, the LPM 
effect does not seriously modify the shower longitudinal profile (in the case of gamma primaries) when the primary 
energy is less than 10^® eV. 

The longitudinal development of the number of muons is also dependent on the LPM effect as shown in figure 

These plots also show the results of simulations made with the LPM effect switched on, but without considering 
the dielectric suppression. The impact of this last effect is of course less important but not small enough to justify not 
considering it, and in fact it is more important in relative terms than in the case of the longitudinal development of 
charged particles (figure ^). It is interesting to describe how the number of muons can be altered when the dielectric 
suppression of electron bremsstrahlung photons is taken into account: In the case of electrons of energy slightly 
below the primary energy, that is, about 10^® eV, the dielectric suppression diminishes the probability of emission 
of photons with energies below ID GeV {y « 10~^), producing a relative enlargement of the number of events with 
photon energies slightly above this limit. Such photons can undergo photonuclear reactions creating pions which in 
turn may decay into muons; and it is clear that if their number is enlarged so will be the number of secondary muons 
(compare the full and dotted lines of figure [lo|b) . 

If the primary particle is not electromagnetic, the influence of the LPM effect will not depend directly on the 
primary energy, but on the energy of the particles that initiate the electromagnetic shower. In the case of showers 
with hadronic primaries, protons for example, the electromagnetic shower is typically started by the gamma rays 
product of decays of 7r° mesons produced after inelastic hadronic processes that take place when the primary collides 
against an atmospheric nucleus. For primary energies larger than 10^^ eV, an inelastic collision hadron-nucleus 
produces hundreds of secondaries, most of them having energies which are f to 4 orders of magnitude smaller than 
the primary energy. Roughly speaking, this means that the electromagnetic shower maximum energy is 2 or 3 orders 
of magnitude below the energy of the primary, and accordingly with the data presented in figure ||, this means 
that the LPM effect will not produce serious distortion in the shower development unless the energy of the primary 
hadron is well above 10^^ eV. This agrees with our results presented in figure Il| which correspond to the longitudinal 
development of 3 x I0^° eV proton showers: The curves corresponding to both LPM "on" and "off" cases do not 
present significant differences. 

Another remarkable feature is the increase of the fluctuations of ATmax and Nnia.x- In figure JlJ (y_3h X^ax (A'max) 
and its fluctuations are plotted as functions of the primary energy, for showers in the range ID ~- 10^^ eV. It shows 
up clearly that: (i) The well-known proportionalities between X^^ax and logiVmax with the logarithm of the primary 
energy [|f],^ do not hold exactly for electromagnetic showers of ultra-high energy when the LPM effect is taken 
into account. However, the mentioned linear relations are valid for proton showers in the whole range of energies 
considered, (ii) As mentioned, the fluctuations of A^max and especially those of A"max become larger as long as the 
primary energy is increased, in comparison with the non-LPM counterparts. The results presented in this section 
correspond to particular cases which are representative of the behavior that should be expected in other general cases. 
Notice that the characteristics of electron initiated showers are very similar to those corresponding to gamma showers 
and the same remark can be done in the case of showers initiated by nuclei in comparison with proton showers. For 
that reason we are not including here any related plots. 



V. CONCLUSIONS 



We have analyzed exhaustively the LPM effect and the dielectric suppression from the theoretical point of view and 
have discussed the practical implementation of the corresponding algorithms for air shower simulations. 

We have studied the different approaches to the LPM effect and conclude that the final results of the Migdal theory 
are best adapted for the numerical treatment while there are no important differences with the results coming from 
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other approaches. 

By means of numerical simulations using AIRES H, we have studied the influence of the LPM effect, including 
dielectric suppression, on the longitudinal development of air showers initiated by ultra-high energy astroparticles. 
As mentioned previously, AIRES is capable of calculating many air shower observables in a realistic environment. 

We have analyzed the influence of both suppression mechanisms in a wide range of primary energies (up to 3 x 10^'' 
eV), and for various primary particles, namely, photons, electrons and protons. 

The main purpose of our work has been to study the influence of the LPM effect on the behavior of electromagnetic 
air showers, for primaries that interact after reaching the Earth's atmosphere. 

Clearly, in the special case of cosmic gamma rays, it would be necessary to take into account the interaction of the 
primary with the geomagnetic field that could take place before entering the atmosphere [p6|-^ . The cosmic photon 
can create electron-positron pairs, which in turn may originate additional magnetic bremsstrahlung photons, and so 
on. As a result, a pre-shower takes place before the particles reach the atmosphere, and consequently the original 
energy is shared among various electromagnetic particles. If the number of particles in the pre-shower is large, the 
influence of the LPM effect will be substantially reduced, in comparison with the case of a single electromagnetic par- 
ticle. Nevertheless, the conversion probability strongly depends on the initial conditions (geomagnetic fleld strength, 
angle between the velocity of the particle and the magnetic fleld, etc.) p^ ], and can be relatively small in certain 
circumstances. Consequently, the pre-shower may eventually not exist or be very small, allowing ultra-high energy 
electromagnetic particles to reach the atmosphere, producing showers like those included in our study. 

We have found that the LPM effect introduces significant modifications on the development of gamma and electron 
air showers if the primary energies are larger than 10^^ eV. The most evident signature of the effect is the shift in the 
position of the maximum of the shower, which moves deeper into the atmosphere with increasing fluctuations when 
the primary energy is enlarged. Our conclusion is that in such cases the effect must be always taken into account in 
realistic simulations of ultra-high energy electromagnetic air showers. It is also important to remark that: (i) The 
influence of the dielectric suppression is not as important as the LPM effect, but large enough to justify having it 
taken into account in any realistic simulation, (ii) For showers with large zenith angles, the suppression that delays 
the shower growth can be not as large as in the case of vertical showers, as explained in section fVl 

We have not found any important effect in proton showers with primary energies up to 10^^ eV. The reason for this, 
as explained in the previous section, is that the electromagnetic shower, where the LPM effect takes place, begins 
later, when the initial energy is shared among the secondary particles and the average energy per particle is then 2-3 
orders of magnitude less than the primary energy. Clearly the same reasoning is valid for nuclei primary cosmic rays. 
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FIG. 1. Functions G(s) (solid line) and $(s) (dashed line), appearing in Migdal theory. The plots correspond to exact 
calculations made by numerical evaluation of equations (^3|) and (p^. 
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FIG. 2. Bremsstrahlung probability according to Migdal theory (kdPjdk). The electron energies are E = 10^* eV (solid 
Une), E = 10^^ eV (dashed line) and E = 10^° eV (dotted line). The medium is air at the vertical depth of Xy = 1000 g/cm^ 
(p = 1.19 kg/m^) (a), and = 50 g/cm^ (p = 78 g/m^) (b). 
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FIG. 3. Bremsstrahlung probability according to Migdal theory (kdP/dk). The probabilities taking (solid line) and not 
taking (dashed line) into account the influence of dielectric suppression are plotted versus logjQ y. The energy of the electron 
is 10^* eV and the medium is air at a vertical depth Xy = 1000 g/cm^ (p = 1.19 kg/m^). 
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FIG. 4. Pair production probability, P(k,u), according to Migdal theory plotted versus u for different values of the photon 
energy: 10^" eV (solid line), 10^^ eV (dashed line) and lO'^^ eV (dotted line). The medium is air at the vertical depth of 
X„ = 1000 g/cm^ (p = 1.19 kg/m^) (a), and X„ = 50 g/cra (p = 78 g/m^) (b). 
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FIG. 5. Comparison between Migdal theory and Blankenbecler and Drell formulation, (a) Absolute probabilities versus y. 
The electron energy is 10^* eV, and the solid (dashed) line represents Migdal theory (Blankenbecler and Drell formulation) . The 
dotted line corresponds to the Blankenbecler and Drell formulation without adding the phase-amplitude correlation correction 
(p9[). (b) Percent relative differences between both formulations, for electron energies 10^* eV (dashed line) and 10^° eV (solid 
line). In all cases the medium is air and the vertical depth is = 1000 g/cm^ (p = 1.19 kg/m'^) (lO'^* eV) and = 100 
g/cm^ (p = 0.156 kg/m^) (10^" eV). 
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FIG. 6. £lpm versus vertical atmospheric depth. 
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FIG. 7. j/diei versus vertical atmospheric depth. 
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FIG. 8. Average fraction of air shower particles with energy larger than -E'lpm/4, (/>), plotted as a function of the vertical 
atmospheric depth. The lines are informal fits to simulation data obtained using AIRES. The solid (dashed) line correspond to 
showers initiated by gamma primaries propagated taking (not taking) into account the LPM effect. In both cases the primary 
energy is 3 x 10^*^ eV and the shower axis is vertical. 




FIG. 9. Longitudinal development of charged particles, obtained by means of Monte Carlo simulations using AIRES, plotted 
versus the vertical atmospheric depth. The primary particle is a photon injected at the top of the atmosphere (100 km above 
sea level) with a zenith angle of 60 degrees. The primary energies arc 10^" cV (a) and 3 x 10^° eV (b). The solid (dashed) 
lines correspond to calculations that include (do not include) the LPM effect. The dotted lines correspond to the LPM case 
but without the dielectric suppression. In some cases the error bars (the RMS errors of the means) were not plotted for clarity; 
they are, in general, smaller or comparable to the represented ones. 
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FIG. 10. Same as figure ^, but for the longitudinal development of the number of muons. 
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FIG. 11. Longitudinal development of charged particles (a) and muons (b) for showers initiated by proton primaries. The 
shower parameters are as in figure bl 
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FIG. 12. Computer simulation results for the shower maximum, Xmax, measured along the shower axis (a), and its fluctua- 
tions (RMS error of the mean) (b) , plotted as function of the primary energy. The solid (dashed) lines correspond to gamma 
primaries taking (not taking) into account the LPM effect. The dotted lines correspond to proton primaries. 
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FIG. 13. Number of charged particles at the shower maximum, A^max, (a) and its fluctuations (RMS error of the mean) (b), 
plotted versus the primary energy. The primaries are gammas in the same conditions as in figure |l^ . 
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TABLE I. The coefficients for the rational expansions of functions G and <& (equations (BSl) and (p4)). 
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